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Abstract Global Circulation Models (GCMs) provide projections for future climate 
warming using a wide variety of highly sophisticated anthropogenic CO2 emissions 
scenarios as input, each based on the evolution of four emissions "drivers": popu- 
lation p, standard of living g, energy productivity (or efficiency) / and energy car- 
bonization c \2\. The range of scenarios considered is extremely broad, however, and 
this is a primary source of forecast uncertainty 1301 . Here, it is shown both theoreti- 
cally and observationally how the evolution of the human system can be considered 
from a surprisingly simple thermodynamic perspective in which it is unnecessary 
to explicitly model two of the emissions drivers: population and standard of living. 
Specifically, the human system grows through a self-perpetuating feedback loop in 
which the consumption rate of primary energy resources stays tied to the historical 
accumulation of global economic production -01 pxg- through a time-independent 
factor of 9.7±0.3 milliwatts per inflation-adjusted 1990 US dollar This important 
constraint, and the fact that / and c have historically varied rather slowly, points 
towards substantially narrowed visions of future emissions scenarios for implemen- 
tation in GCMs. 

Keywords carbon dioxide emissions • growth model ■ evolution • thermodynamics • 
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1 Introduction 

GCM projections of 21st century climate change use prognostic trajectories for car- 
bon dioxide (CO2) emission fluxes developed by the International Panel on Climate 
Change (IPCC) Special Report on Emissions Scenarios (SRES) |2|. These provide a 
range of timelines, each designed to show how a given set of decisions might corre- 
spond to a particular atmospheric CO2 trajectory. SRES models are highly sophisti- 
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cated, and contain numerous interactive components, each designed to reflect a real- 
istic range of societal dynamic behavior. 

For tractability, IPCC SRES models express the primary drivers of growth in CO2 
emissions E in terms of human population p, primary energy consumption a and 
real (or inflation-adjusted) economic production P through E=pxgxixc where 
g = P/p represents the real economic production per person and /= l/f — a/P 
represents the "energy intensity" of real economic production, or alternatively, the 
inverse of its "energy productivity" /, and c ~ E /a the carbonization of the energy 
supply [,19J . Expressed in a prognostic form, emissions grow at a rate given by |fT9l 

m 

dlnE dlnp ding dlnf dine 
dt dt dt dt dt 

Differences among SRES emissions trajectories depend on how society is as- 
sumed to manage such issues as population control, energy efficiency, and a switch 
to non-C02 emitting energy resources. Currently, the range of possible futures con- 
sidered is extremely broad. In fact, uncertainty in the degree of surface warming over 
the next century is determined as much by the range of SRES scenarios as by climate 
physics itself ||30l . 

In this paper I propose that by using a straight-forward thermodynamic approach 
it may be possible to substantially constrain plausible timelines for future anthro- 
pogenic CO2 emission rates. 



2 A thermodynamic growth model 

2.1 A heat engine 

The starting point is recognition that general thermodynamic laws require that all 
systems, even those that are living, evolve through a spontaneous conversion of envi- 
ronmental potential energy into some less available form, often termed "heat" ll27l [8l 

mm- 

Specifically, consider a system drawn in Fig. [T]consisting of some entity and its 
environment, separated by some arbitrarily defined permeable interface at a fixed 
temperature and pressure p (i.e., at constant energy density or an isentrope). As 
a whole, the system is in contact (through radiation, convection or conduction) with 
colder, lower energy density surroundings at temperature T < T^. The interface be- 
tween the entity and the environment represents a "step", with its height represented 
by a Gibbs energy potential AG. Available potential energy in the environment is 
converted at rate a — a AG into some less available form through the transfer of mat- 
ter across the interface. The system-specific constant coefficient a is an intensive 
quantity that defines the particular physics of "availability" for the systemQ. 

Effectively, the system operates as a form of "heat engine". The familiar textbook 
heat engine has the engine consume energy at rate a to do "work" at rate w to con- 
tribute to the potential of some outside agency while releasing waste "heat" at rate 

' To take an electrostatic analogue, zi G is a voltage difference, a a current, and 1 /a is the resistance, in 
which case the relevant physics defining "availability" of energy is the material's conductivity. 
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Fig. 1 Illustration of an evolving system bounding some entity and its environment, as separated by a 
permeable interface at constant temperature and pressure. The interface maintains an energy potential 
AG{Ts,p) so that the system as a whole is able to convert available energy at rate a into work w with 
efficiency e = wla, and "heat" at rate a — w. Heat is voided to the system's colder surroundings; work 
grows the interface at rate w = d(AG)/dt. Because the interface potential is related to energy consumption 
through a = aAG, where a is an engine specific constant coefficient, what is defined is a positive feedback 
loop in which, through work, a and ziG evolve logarithmically at rate dlna/dt = dla{AG) /dt = r] = ea. 
Here 77 can be considered a feedback efficiency or rate of return. 



a~w |f35l . While the definition of work is clear for an industrial steam engine raising 
the gravitational potential of a steel beam, for example, the choice of what qualifies 
as work is quite subjective. In fact, all energy transfers a act to increase the potential 
of something. Work is simply the raised potential of interest. Heat is the remainder. 

With reference to Fig. [T] work is subjectively defined with respect to the inter- 
nal energy potential AG of the interface separating the environment and the entity. 
Thus, through consumption of available energy at rate a, the value of AG evolves 
at rate w — dAG/dt with heat engine efficiency e = w/a. Meanwhile, heat is lost 
spontaneously to the colder surroundings at rate a — w. The Second Law of Ther- 
modynamics requires that heat production a — w > 0, in which case e must be less 
than unity. Thus, the existence of a potential difference A G between the entity and 
its environment entails decay of the available potential of the universe as a whole (or 
equivalently, an increase in its entropy) 1.35] . 

The advantage of the above thermodynamic setup is that it allows for spontaneous 
evolution of the entity, and, as will be shown, it can be applied more specifically 
to the evolution of civilization. Because work is internal, a feedback loop causes 
the interface to exponentially grow or decay: the existence of AG requires energy 
consumption at rate a = aAG; in turn, this corresponds to work being done at rate 
w = ea, which then adds to the internal potential at rate dAG/dt = w, closing the 
loop. If work increases the magnitude of AG, then the interface separating the entity 
and its environment bootstraps itself to a higher level. Then the system as a whole 
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evolves to higher levels of energy consumption a through 

da d(AG) 

— = a— -^aw~aea = r\a (2) 

dt dt 

where r\ is effectively a "rate of return" representing the efficiency of the feedback 
on energy consumption a. Note that, perhaps counter- intuitively, higher energy effi- 
ciency e corresponds to higher values of Tj, and therefore more rapidly exponential 
evolution of energy consumption a and heat production a~w. 

In AppendixlAl the nature of the feedback efficiency 7] is defined more precisely. 
It is shown that the interface AG can be separated into h material units, each asso- 
ciated with the same potential energy at fixed temperature and pressure of Apt. AG 
results in a flux of material across the interface at rate dn/dt ^ a/ Ajj.. If the net flow 
is from the environment to the entity, a portion of material that diffuses across the in- 
terface at rate dn/dt then contributes to interface growth at rate dh/dt. The feedback 
efficiency 77 is the logarithmic form for this material rate of growth 

d\nh 

A concrete example that might be particularly easy to relate to is the growth of 
a young child. As an entity, the child consumes the accessible energy contained in 
food from the environment in proportion to some measure of the child's size. This 
rate of consumption a ~ a AG - perhaps about 50 Watts - enables the child to do 
"work" w = dAG/dt with energy efficiency e = w/a, incorporating the water and 
nutrients contained in food into its structure in order to extend the material interface 
n separating it from its environment. The child maintains homeostasis because "heat" 
can eventually radiate to space at rate a — w at a relatively cold planetary blackbody 
temperature of about 255 K. Material waste is also produced once the useful chemical 
potential of the nutrition has been extracted, for example as carbohydrates are con- 
verted to exhaled CO2. Through a feedback loop, if w > 0, the child and its energy 
consumption grow logarithmically at a rate rj — dlna/dt. Of course, in an energy 
poor environment there might not be sufficient nutrition, in which case w < and the 
feedback efficiency rj is negative. But, assuming the child reaches adulthood, growth 
tends towards a balance between energy consumption and heat production, and 77 
tends to zero. 



3 Analog for the growth of civiUzation and its emissions 

The argument now is that the thermodynamic growth model described above, just 
as it can be applied to a child's growth, can also be extended to the human system 
in its entirety, as defined by civilization and its known environmental reservoirs. As 
with the child, an interface potential A G between civilization and its primary energy 
resources enables energy to be consumed at rate a. This allows work to be done with 
efficiency e and at rate w to grow the interface potential AG through incorporation 
of environmental matter (e.g., biomass and minerals). Simultaneously, through con- 
vection and radiation, heat is lost to space at rate a — w at the planetary blackbody 
temperature. Also, unused material waste accumulates in the environment. 
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Of course, for civilization, "food" includes the chemical and nuclear bonds in 
oil, coal and uranium, combined with mineral matter from the Earth's crust. These 
build structures that include not just human bones, flesh, blood and nerves, but cities, 
roads, shipping and telecommunications. However, insofar as the thermodynamics is 
concerned, the difference between the child and civilization is really only a matter of 
complexity and scale. In either case, as part of a single energy-consuming organism, 
all organism elements contribute to an interface with environmental reservoirs that 
enables net available energy transfer to the organism at rate a. 

While a precise definition of civilization is arbitrary, civilization is most com- 
monly quantified in purely fiscal terms. Thus, the goal here is to examine whether it 
is possible to link fiscal quantities to the more thermodynamic model defined above. 
To this end, an argument can be made that, if what physically distinguishes civiliza- 
tion from its environment is some thermodynamic potential AG at constant tempera- 
ture and pressure, civilization implicitly assigns inflation-adjusted (or real) monetary 
value to what AG enables - the total rate of energy consumption a. To borrow a 
phrase, "money is power" because, if all current exothermic processes supporting 
civilization were to suddenly cease such that a equalled zero, all civilization would 
become worthless; it would no longer be associated with a non-equilibrium level of 
potential energy AG^ a/a. Simply, there would be no definable material interface n 
between civilization and its environment. 

As an example, the potential energy in oil combustion is valuable, but only to 
the extent it that it can interact with the interface separating civilization from its en- 
vironment. It has zero value if it burns wastefully in the desert, and zero value in 
its unavailable chemical and nuclear bonds. From society's perspective, any societal 
element, whether living or synthetic, only has value to the extent it is able to op- 
erate in synergy with all other elements to define an interface with environmental 
available energy. An unavailable road from nowhere to nowhere is just pavement on 
the ground. But the same road between two cities is part of a larger organism that 
works collectively at net rate w to grow access to the primary energy supplies that 
civilization requires. 

The mathematical expression of the above argument is that global primary energy 
consumption a is related to global value C through a constant factor X 



Thus, the economic representation of the evolving heat engine given by Eq.|2]is 



a = XC 



(4) 



dC 1 da a d{AG) a 




(5) 



dt X dt X dt X 



—w 



or in purely economic terms 




^c 



(6) 



where, C (units value) grows through the real (inflation-adjusted) economic produc- 
tion rate P (units currency per year). The thermodynamic feedback efficiency Tj is an 
economic rate of return on C. 
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Eq.|5] implies that real economic production P is, perhaps rather intuitively, only 
a measure of thermodynamic work w through coefficient a/ A. Expressed in integral 
form 

C{t)= f P (f ') dt' - ^ /' w (f ') dt' (7) 

70 A Jo 

Thus, this growth model is a statement that the rate of return (or feedback efficiency) 
r\ — aeon economic value C is a consequence of doing thermodynamic work w with 
efficiency £ to grow the interface between civilization and environmental resources. 
By growing the interface, civilization is able to draw more energy a, and do more 
work w, thereby closing the loop. 

Now, returning to frameworks for CO2 emissions forecasts, the SRES definition 
for energy productivity f ~ P/a can be understood in light of the above. If Eq.|2]is 
combined with Eqs.|5]and|6] this yields the basic relations 

r\=P/C = ae = Xf (8) 

Therefore, energy productivity f — P/a\s related to the heat engine thermodynamic 
efficiency e — w/a and a fiscal expression for the feedback efficiency t] — P/C 
through the fixed, intrinsic quantities a and A . 



- 

dt 

V 

A' 

t + 

da 

dt < — 

V 

Fig. 2 Schematic extending Fig.[T]to relate energy consumption a by the human system to economic value 
C and CO2 emissions E. Black arrows point in the direction of the product, red arrows in the direction of 
the integral over time. Work is done at rate w to enable energy consumption a to grow at rate da/dt = r)a, 
where r) is the feedback efficiency of a heat engine representing the system. The economy has a fixed 
relationship to energy consumption thi'ough a = AC, where C is civilization's historical accumulation of 
real (inflation-adjusted) economic production of economic value P = dC/dt (units currency), and A is 
an intrinsic constant of proportionality. Thus, CO2 emissions are related to economic production through 
E = Xc Jq P (t')dt' , where c is the carbon content of energy in the fuel supply. 
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A schematic illustrating the economic growth model is shown in Fig. |2] A dis- 
cussion of how it relates to more orthodox economic approaches is contained in Ap- 
pendix [B] where it is shown how such traditional economic concepts as inflation, 
savings and capital depreciation can be interpreted within a thermodynamic context. 
A straightforward consequence of Eqs.|5]and|6]is that the rate of growth of the global 
economy obeys the simple relation 

dlnP dlnri 

Interestingly, the approach is of identical mathematical form to one often used to 
successfully model growth of vegetation, where vegetative "value" C refers not to 
money but instead to biomass, and P to the net primary productivity 031I I18I. Presum- 
ably, biological organisms must also maintain a high potential interface with respect 
to their environment, enabling them to consume energy, produce heat and waste, and 
do work to incorporate the matter that enables them to grow |[33l . Thermodynamic 
laws are fully general. 

A difference between plants and civilization is that plant waste includes CO2 
that is recyclable, whereas the global economy creates most CO2 from fossil-carbon, 
much of which accumulates in the atmosphere. From Fig. |2l CO2 emissions can be 
represented simply through 



E{t)^XcC = Xc j'^P{t')dt' 



(10) 



Present-day emissions are determined by past accumulation of real economic pro- 
duction and the current carbonization of the energy supply. Current emissions growth 
rates are given by 

d\\\E d\nc 

Eq.[TT]is more simple and physical than the expression for drivers in SRES forecasts 
given by Eq.[T] 



4 Evaluation 

The preceding discussion argues for a direct theoretical link between anthropogenic 
emissions and basic thermodynamics. But is the argument observationally supported? 
The expression for emissions growth, Eq. [TT] rests on the premise that there exists 
an intrinsic quantity A representing how the historical accumulation of economic 
production C is supported by a rate of energy consumption a (Eq. [Hi. If A is not 
constant with time, then the thermodynamic framework is false. 

I examine this proposition now using statistics for the combination of world en- 
ergy production a |1| and real global economic production P |3| (expressed here 
in fixed 1990 US dollars) for the 36 year interval between 1970 to 2005 for which 
these statistics are currently available. The time series for accumulated global eco- 
nomic value C = jQP{t')dt' is estimated using sporadic calculations of P that have 
been ascertained for select years over the past two millennia |fT6l in combination with 
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Fig. 3 Estimates of gross world product P in market exchange rate, 1990 US dollars and economic value C, 
defined by P = dC/dt. Also shown are recent global primary energy consumption a, the ratio X=a/C, and 
the feedback efficiency 7] = P/C. Dashed lines correspond to extrapolations based on assuming X = 9.7 
mWper 1990 US dollar. 



more recent annual records 1 3 1 to create a two-millennia yearly time-series in P [See 
AppendixICl. Estimates of P and C and their ratio rj = P/C are shown in Fig.|3] 

Figure [3] shows that, over a period between 1970 and 2005, the ratio X (f ) ^ a/C 
maintained a nearly constant value of 0.306 exajoules per trillion 1990 US dollars per 
year, or alternatively 9.7 milliwatts per 1990 dollar. Corrected for autocorrelation in 
the time-series, the observational uncertainty at the 95% confidence level is just ±0.3 
milliwatts per 1990 dollar. The simplest interpretation is that this result supports the 
cornerstone hypothesis given by Eq. lU the historical accumulation of real economic 
value through real economic production is maintained by continuous primary energy 
consumption; the relationship between value and rates of energy consumption is a 
constant parameter. 

Of course it is possible that this observed result only holds over the 36-year pe- 
riod for which global energy consumption statistics are available, but it is expected 
theoretically; the period examined covers over half of total historical growth in a and 
C, and two thirds of P; and, the observational uncertainty is small enough to plausibly 
reflect errors or noise in historical data. For example, new primary energy production 
(what has been measured) only reflects new primary energy consumption (what is 
theoretically relevant) in the average, not the instant. 
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5 Drivers of emissions growth 

The existence of a fixed relationship between energy production and accumulated real 
economic production simplifies the number of drivers required for CO2 emissions 
forecasts. To see how, the SRES emissions growth equation Eq. [T]can be equated 
with the more thermodynamic expression given by Eq. [TT] Since both expressions 
rely on exogenous expressions for carbonization growth dlnc/dt, this is effectively 
a comparison of expressions for growth in energy consumption dlna/dt 



^ + ^-^ = r,^A/ (12) 
dt dt dt 



SRES models consider population p and standard of living g and energy productivity 
/ as the key "drivers" of energy consumption growth, but the "driver" concept can be 
misleading when, at a very basic level, feedback determines how p, g, and / are inter- 
related. Eq. [12] demonstrates that growth in p and g is fundamentally constrained by 
the sum of the current state of the feedback efficiency 77 = A/ and its rate of change 
dlnf/dt. Therefore, knowledge of the behavior of only one parameter, /, is required 
for forecasts of energy consumption growth, rather than each of /, p and g. 

So, perhaps surprisingly, changes in population and standard of living might best 
be considered as only a response to energy efficiency. As part of a heat engine, cre- 
ating people and their lifestyles requires energy consumption. Doing so efficiently 
merely serves to bootstrap civilization into a more consumptive (and productive) 
state. Likely, society has traditionally praised energy efficiency gains for precisely 
this reason. As summarized by Eq. |2] energy efficiency gains accelerate rather than 
slow energy consumption II12II29I , contrary to what is commonly assumed (20] 

This does not mean that consideration of population and standard of living in 
SRES models is invalid, of course, only that their evolution must be consistent with 
thermodynamic constraints (Eq. \l7\. Fig. |4] shows a time-series comparing 15-year 
hind-casts with observations, evaluated for both sides of Eq. [12] in each of the years 
1980 to 1990. Hind-casts based on growth of p, g and / apply straight-forward per- 
sistence in trends from the prior 10 years, i.e., dlna/dt ~ dlnp/dt + dlng/dt — 
dlnf/dt. For comparison, hind-casts based on the zeroth-order thermodynamic ex- 
pression for economic growth need only employ evaluations of the current-year state 
of Tj = P/C. 

Since both approaches reflect already realized thermodynamic constraints, both 
"persistence" and a thermodynamic model give hind-casts that reproduce observed 
trends with comparable accuracy. Notably, the thermodynamic model provides a 
hind-cast for average growth between 1990 and 2005 that is within just 0.1 %/yr 
of observed growth rates. A prior study found that this level of accuracy was only at- 
tained by a particular "worst-case" SRES model for this particular time period ll24ll . 
What is important here is that, for the purpose of future forecasts, the thermodynamic 
approach is accurate while being both simpler and more physical than using persis- 
tence or sophisticated SRES models. 
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Fig. 4 Difference between predicted and actual growth rates in global energy consumption a, based on 
predictions derived using persistence in trends (red) and the zeroth-order model presented here (green). 
Trend persistence is calculated using the prior ten years as a basis; the zeroth-order model is based only on 
current year calculations of the feedback efficiency r\ . Hind-casts are compared with actual observed rates 
for a period covering the following 15 years. 



6 Considerations for modeling future emissions scenarios 



An advantage of appealing to energy efficiency in forecasts of CO2 emissions is that 
77 = A/ tends to vary rather slowly. Since 1970, growth of T] = P/C has climbed 
from 1.4 % per year to 2.1 % per year in 2005 (Fig.O, corresponding to an e-folding 
time-scale Tjj = \/{d\nr\/dt) of approximately 100 years. Expressed in terms of 
time-series analysis, Tj is highly "reddened", because it is an integrator of drj/dt. 
Alternatively, and in more fiscal terms, global economic value C (and hence energy 
consumption a = AC) varies slowly because it is an integrator of economic production 
P = dC/dt (Fig. [3]). The present and future are influenced by even the most distant 
past, and the past cannot be erased. 

The carbonization of the energy supply c, is changing even more slowly (see 
Appendix iDb with a time scale of about 300 years. What this means is that future 
emissions rates for CO2 are most strongly influenced by the current state of T] . As 
a zeroth-order assumption, it is reasonable to assume persistence in 7], meaning that 
over time-scales much less than T-q , future emissions are unlikely to depart substan- 
tially from the recent growth rate of 2. 1 % per year. 

More accurate forecasts of energy consumption and CO2 emissions rates will 
require an understanding of how t] itself evolves. Assuming c is a constant, positive 
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values of dlnrj /dt imply super-exponential growtl0of CO? emissions E ||2TI . The 
solution for Eq. (TT] starting at some time f,-, is 



E = Ej exp 



(13) 



Note, that growth condenses to the single exponential form in the limit of Zif <C T,j . 

However, in the long-term, assuming persistence in even dlnrj /dt is an over- 
simplification since T-q itself evolves. History shows bursts in efficiency growth, no- 
tably around 1880 and 1950, perhaps when important new discoveries of energy 
reservoirs made the past less relevant (Fig. O. But, in both cases, the initial burst 
in Tj eventually tapered. After 1950, the time-scale Tri, changed from just 30 years 
between 1950 and 1970, to 67 years between 1970 and 1990, and 120 years between 
1990 and 2005. Plausibly, 1/t,j = dlnrj /dt will eventually cross zero and turn nega- 
tive, implying sub-exponential growth in emissions E (Eq.[T3Tl. 

Unfortunately, if —li In 77 /dt is ever greater than rj over the long term, while emis- 
sions growth may be significantly slowed, what is implied is a real global economy 
that is shrinking (Eq. |9]l. Robust multi-decadal forecasts of emissions E, and its re- 
lationship to economic production P, require a first principles thermodynamic model 
for how TJ changes with time. Assuming 1880 and 1950 were indeed associated with 
discovery of new energy reservoirs, this would suggest the problem is fundamentally 
geological, and that higher-order moments of rj reflect rates of reservoir discovery 
and depletion. 

Understood thermodynamically, the transfer of energy at rate a across the inter- 
face between energy reservoirs and civilization reflects a balance. On one hand, the 
transfer grows civilization, and increases the physical size of the interface AG. At the 
same time, however, it depletes the reservoirs, and this decreases AG. The sign and 
magnitude of the rate of work w — dAG/dt, and therefore 77, depends on the relative 
strengths of these two forces. 



7 Mitigation 

The premise behind mitigation is that there are "drivers" of emissions rates that can be 
meaningfully controlled through policy. As shown above, it is not clear that the driver 
concept is in fact meaningful. Rather, it appears that drivers in the Kaya Identity are 
merely a thermodynamic response to the current value of rj . 

At least this appears to be true for population p and standard of living g. It is 
not yet clear whether it applies to the current carbonization of the economy c. An 
interesting result that can be derived from Eq. [TO] using the values for X in Fig. [3] 
and c in Fi g . 15) (Appendix IdTi is that the "carbon footprint" of civilization in recent 
decades reflects a simple relationship between the rate of global carbon emissions 
E and the accumulation over history of real global value C. The coefficient is Ac = 
5.2±0.2 MtC per year, per trillion 1990 US dollars of global economic value. 

- Similar super-exponential growth behavior has been obser'ved previously at a more local level, in the 
characteristics of cities 1 6 1 . 
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To take the result further, Eq. [TT]points towards a non-dimensional number 

representing the relationship between the global economy's rate of de-carbonization, 
—dlnc/dt, and its rate of return, rj = A/. If 5 > 1, dlnE /dt < 0, and emissions are 
stabilized or declining. 

To reach stabilization, what is required is decarbonization that is at least as fast 
as the economy's rate of return. Taking the 2005 value for rj of 2. 1 % per year, stabi- 
hzation of emissions would require an equivalent or greater rate of decarbonization. 
2.1% of current annual energy production corresponds to an annual addition of ap- 
proximately 300 GW of new non-carbon emitting power capacity - approximately 
one new nuclear power plant per day. 



8 Conclusions 

The physics incorporated into GCM representations of the land, oceans and atmo- 
sphere is required to adhere to universal thermodynamic laws. Ideally, the CO2 emis- 
sions models meant for implementation in GCM projections of climate change should 
do so as well. Fortunately, it appears that appealing to thermodynamic principles may 
lead to a substantially constrained range of possible emissions scenarios. If civiliza- 
tion is considered at a global level, it turns out there is no explicit need to consider 
people or their lifestyles in order to forecast future energy consumption. At civiliza- 
tion's core there is a single constant factor, A = 9.7 ± 0.3 mW per inflation-adjusted 
1990 dollar, that ties the global economy to simple physical principles. Viewed from 
this perspective, civilization evolves in a spontaneous feedback loop maintained only 
by energy consumption and incorporation of environmental matter. 

Because the current state of the system, by nature, is tied to its unchangeable past, 
it looks unlikely that there will be any substantial near-term departure from recently 
observed acceleration in CO2 emission rates. For predictions over the longer term, 
however, what is required is thermodynamically based models for how rates of car- 
bonization and energy efficiency evolve. To this end, these rates are almost certainly 
constrained by the size and availability of environmental resource reservoirs. Previ- 
ously, such factors have been shown to be primary constraints in the evolution of 
species |33 ,34|. Extending these principles to civilization, emissions models might 
be simplified further yet. 
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Appendix 

A Material transfer 

To understand the details of the heat engine described in this article in a bit more detail, it is helpful to look 
more explicitly at what constitutes available energy and work for this case. The Gibbs energy potential of 
matter can be expressed as Y.i'^it^i [T-.p), where «, refers to the number of the species (' with specific chem- 
ical potential jli(T,p) I35I13I (the "chemical" potential, rather confusingly, can always be generalized 
where relevant to incorporate the potential exothermic energy in nuclear bonds). The interface separating 
the entity and its environment may be composed of matter in many forms. However, a simplifying argu- 
ment can be made that the potential difference zi G can be split up into h arbitrary units of matter, each unit 
candying an identical available potential of A^{T,,p) 

AG = nA^{T„p) (15) 

Likewise, the net flux of material between the environment and the entity at rate dn/dt requires energy 
consumption by the system as a whole at rate 

a=^An{T„p) (16) 

But, since a = azi G, it also holds that 

a = anAil(Ts,p) (17) 

Combined, Eqs. ll6l and ll7l imply that the intensive quantity a = l/« (dn/dt) is determined by the partic- 
ular physics relating the amount of high potential matter along the interface to the flux of matter across it|3 

The evolution of energy consumption by the system, da/dt, is related to its rate of doing work 
through aw (Eq. |2j, but more specifically to the interface's material growth. Since work is defined by 
w = d (AG) I dt , and the interface temperature and pressure are fixed, it follows that the potential defining 
the interface between the entity and its environment evolves at rate 

w=^An{T„p) (18) 

Work is positive if the material interface grows. Expressed in terms of the rate of energy consumption, 

dit 1 dlnii 1 da 

w=—a = — a = — — (19) 

dn a dt a dt 

Because the heat engine efficiency is given by e = w/a, this leads to the result that 

1 d\nh 

e=dndn = (20) 

a dt 

As a feedback loop, Eqs.l2land l20l can be combined to show that the feedback efficiency 7] is 

d\nn 

Thus, r\ expresses the logarithmic growth of the number of elements defining the interface between the 
entity and its environment. If the net flow is from the environment to the entity, the portion of material that 
diffuses across the interface and does not contribute to interface growth is returned to the environment as 
waste. 



' It is straightforward to show that for the special case of Maxwellian diffusion along a concentration 
gradient to, for example, a cloud droplet or snow flake |23 1, evolution of a and n is determined not by the 
surface area of the interface (as might initially seem more intuitive) but rather by a length dimension. In 
this case, a is determined by the product of the diffusivity of vapor in air and the area density of vapor at 
saturation. 
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B Comparison with traditional economic models 

Economic studies normally separate production into two components: a fraction i representing a sav- 
ings, or investment; and a fraction (1 — .s) representing private and government "consumption". Models 
represent the nominal growth in "capital" K (units currency) as the difference between the portion .v of 
production P (units currency per time) that is a savings or investment, and capital depreciation at rate 7 

— = {P-W)-yK = sP-yK (22) 

where individual and government consumption is represented by W = ( 1 — .v) P. 

In return, according to some functional form, labor L (units worker hours) employs capital K (units 
currency) to generate further production P (units currency per time). For the sake of illustration, a com- 
monly used representation is the Cobb-Douglas production function 

P = AK"L^'°- (23) 

where A, the "total factor productivity", is a compensating factor designed to account for any residual 
unaccounted for by K and L. The exponent a is empirically determined. The Solow Growth Model 1281 
expresses the prognostic form for Ea. l23l as 

dlnP d\nA d\nK , , dlnL 

—— = —— + «— —+l-a—— 24 
dt dt dt dt 

Commonly, the term d\n A j dt is interpreted to represent technological progress. 

There have been criticisms raised of the Solow Model because it makes no explicit reference to natural 
resources 110151 . One suggested remedy has been to incorporate primary energy consumption into Eg. 1231 
as a complement to labor or capital 1 25 26], in which case 

P=(AKK)"(ALLf{A,a)'-''-l' (25) 

where, again, a is energy consumption, a and /3 are empirically determined, and the subscripts for A refer 
to respective technological progress. 

Now, by comparison, in the thermodynamic economic growth model introduced here, real (inflation- 
adjusted) economic production P (units currency per time) and global value C (units currency) are, respec- 
tively, fiscal representations of net thermodynamic work and the rate of consumption of available primary 
energy resources. The economic growth model described by Eqs|5]and|6]is given by the value production 
function 

P=r\C (26) 

where 7] is the feedback efficiency representing a rate of return due to thermodynamic work by the system 
on the system. The equation for growth of value is 

dC/dt = P (11) 

Note that while C is analogous to capital K in Eq.|22] since C = a/A, it is a more explicitly thermodynamic 
expression of value. 

So, to put the above in context of standard economic production functions, Eq.|26]can be considered 
to be a simplification of Eq. [25] The representation of economic value C employed here is a substitution 
of the combination of traditionally defined capital K and labor L in Eq. |25] such that a = I and /3 = 0, 
and An = r\ . Alternatively, since C is itself only a monetary representation of the rate of primary energy 
consumption a through a = AC, it could equally be stated that a = p = 0, and A„ = 7] /A . 

Note that the thermodynamic production function (Eq. l26K unlike more standard formulations (Eq. 
I23t . has the mathematical advantage of being dimensionally self-consistent, as it does not need to appeal 
to non-integer exponents a and j3 of dimensional terms (such as L and K), as fitted to a specific set of 
circumstances, and with no certain application to different economic regimes. 

It might be argued, however, that the model introduced here fails by leaving no room for either con- 
sumption W or depreciation yK, two central components of the standard economic growth equation (Eq. 
I22t . Offhand, this seems reasonable because, certainly, some portion of economic production must be con- 
sumed, at least in order to maintain economic capital against depreciation or decay: buildings crumble; 
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bodies must be maintained; old technology becomes obsolete; as does past acquisition of human skills and 
knowledge. 

But these concerns can be resolved once it is recognized that the equations derived for this study 
are intended to apply only to real, inflation-adjusted production P, and not nominal production P. To 
demonstrate, assume inflation is positive, in which case nominal value C grows faster than real value by 
some fractional rate of real value, 7 

dC dC 

^ = ^ + ^ '''' 
Since it has been argued here that dC/dt = P, this leads to 

dC ~ 

-^-P-rC (29) 

in which case, the source of real value is nominal production, and the corresponding sink for real value 
occurs at rate 7. So, in fact, Eq. [29] illustrates that Eg. 1 27 1 does account for depreciation through the term 
7C, and is thus similar to the depreciation term yK in the standard growth equation for capital (Eq. l22t . 
While depreciation is implicit when the growth equations are expressed in real, inflation-adjusted terms, 
depreciation is explicit when they are expressed in nominal terms. 

In fact, it is interesting to see what the value decay rate 7 represents. Again, because dC/dt = P, this 
means Eg. 1 27 1 leads to the statement P — P = yC. Alternatively, when nominal production is expressed in 
energy consumption co-ordinates through substitution of the expression a = AC 

P-P = yC=ja (30) 

Compare this to an equivalent expression derived for real production P = {r\ I X)a (Eq.|5). The implication 
here is that economic inflation, the difference between nominal and real production, is a consequence of 
the spontaneous decay or depreciation of total economic value C at rate yC. Put another way, since P = r)C, 
the ratio 7/(7+ J] ) is the fraction of nominal production P that, unlike P, does not return itself as a real 
addition to total value C. In thermodynamic terms, value depreciation yC is an energy barrier that must first 
be crossed for real production to occur. If it is not, the perspective of civilization is that nominal production 
may be positive, but real production is negative. Net work vv is done on civilization by the environment, 
rather than the reverse. Whenever this occurs, the interface between civilization and its environment AG 
decays. 

There is also a consumption term in the traditional expression for capital growth Eg. l22l that is not 
present in the thermodynamic expression for total value growth Eq.|27] As it is normally defined, consump- 
tion is the portion of economic production that does not represent an investment or savings in traditional 
representations of capital K. By contrast, in the thermodynamic model, effectively all real production is 
an "investment" in total economic value C. While a portion of nominal production or nominal work may 
merely serve to offset depreciation of C as described above, all of the remainder adds to the total. Real 
production is net production. 

To illustrate, the construction of coal mines and power plants clearly represents an investment in 
economic value in either framework. A less obvious, although functionally equivalent example, is food 
consumption. In standard representations, food would be "consumed" by households and not contribute 
to their value. However, the available chemical potential in food consumption dn/dtAfl (Eg. I16t also 
maintains and improves that household's capacity to fuilher consume energy and do work by supporting 
its internal potential energy AG (Eq.\Q. Of course, the consumption of an ordinary sandwich may only 
offset a body and mind against decay from "heat" loss, maintaining its internal potential such that it can 
continue to consume energy at the same rate it has in the past (in which case the real production rate P and 
net work rate w is zero since P = yC). The added value of a really good, if more expensive, sandwich is its 
capacity to facilitate real production and new energy consumption above and beyond decay (in which case 
real production is greater than zero and P > yC). The addition to total global value C (and internal potential 
A G) may deiive from a heightened sense of personal well-being and an increased desire to productively 
interact with the rest of civilization in order to afford such sandwiches. 

It is worth noting that a primary conclusion of this paper, that feedback loops in the economic system 
mean that increases in energy efficiency correspond to greater energy consumption, has been reached 
previously by some economists, albeit in a less explicitly physical form than presented here 1 5 25 26 . 14l 
[7j 4 22 9 1 1 1. Although the concept was first introduced by W. Stanley Jevons over a century ago 1 12], the 
extent of energy efficiency "rebound" or "backfire" remains disputed |29|, with no consensus on how it 
should be quantified on the global scales relevant to forecasts of climate change from anthropogenic CO2 . 
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C Materials and methods for time series estimates 

us Department of Energy statistics for global primary energy production Q] include fossil fuel, hydro- 
electric, nuclear, geothermal, wind, solar, and biomass sources. It is assumed here that production and 
consumption rates are, at least on average, equivalent. United Nations time series for world economic pro- 
duction 1 3 | represent the total gross domestic product of all countries, adjusted for inflation and market 
exchange rates to fixed 1990 US dollars. Statistics for CO2 emissions ai'e obtained from the Carbon Diox- 
ide Information Analysis Center 1 17 1. Rather than looking at nations or sectors, only global quantities are 
considered here because, at this level, atmospheric CO2 is well-mixed, and international markets make 
details in economic trade unimportant. 

Gross World Product estimates in 1990 market exchange rate dollars ai'e available for each year since 
1970 |3 1. Long-term but intermittent historical estimates are available for the years 1 to 1992 CE 1 16|. The 
latter data set is expressed in Geary Khamis purchasing power parity (PPP) 1990 US dollars. In general, 
the motivation for expressing valuation in PPP instead of exchange rate dollars is to account for disparities 
in product valuation that exist between countries. In PPP dollars, product valuation is equalized according 
to its apparent contribution to standard of living. Countries with a low standard of living tend to have a 
relatively high gross domestic product when expressed in PPP rather than market exchange rate dollars 
because equivalent products and services tend to be less expensive. 

However, because the focus of this study is energy production and associated CO2 emissions, rather 
than national standard of living, it is historical records of market exchange rate valuations that are used. 
Exchange rate measures of production P are assumed to most accurately reflect the total energy costs 
associated with manifesting products and services in the respective nations where they are consumed. 

To account for any discrepancy between PPP and exchange rate estimates in historical records for 
economic production P, market exchange rate data from 1970 onwards is used to devise a time-dependent 
correction factor Tl to be applied to PPP records such that 7t = PPP/exchangerate. For the period 1970 to 
1992, during which both PPP and market exchange rate estimates of P are available, the fitted value for 
;r is ;r = 1 + 0.258 exp [{f — 1998) /73]. This correction factor can be extrapolated and applied to all PPP 
data between the years and 1969. For the period from 1970 onwards, measured exchange rate values 
are used. Because the historical estimates of P in PPP dollars are increasingly sparse with distance back 
in time (e.g. there are only three data points for the period 1 to 1500 CE), the corrected dataset for P is 
mapped to a yearly distribution using a cubic spline fit. 

Estimates of economic value C represent an accumulation of economic production P over time since 
1 CE, i.e., C{t) = C(l) + S[P{t') dt' . To estimate a value for C(l), it is assumed that the ratio of popu- 
lation to economic value in 1 CE. is equivalent to the average value between 1 CE and the threshold of 
the industrial revolution circa 1700 CE. From historical population statistics 1 16], the associated iterative 
solution for C(l) is 120 trillion 1990 U.S. dollars. For comparison, the estimated value of C in 2005 CE 
is 1580 trillion 1990 US dollars (Fig. 3). Although, off-hand, this value for C(l) seems surprisingly high, 
it is still very small compared to current day values, so the derived value of A presented in this paper is 
relatively insensitive to eiTors in its estimate. 



D Summary of observed growtli rates between 1970 and 2005 

A summary of observed growth rates in global world real production P, carbon dioxide emission rates E, 
feedback efficiency r\ and carbon dioxide emission intensity c between 1970 and 2005 is provided in Fig. 
|5]and in Table[T] along with relevant equations based on the thermodynamic model described here. 
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Fig. 5 For the period 1970 to 2005, trajectories in real global world production P and carbon dioxide 
emissions E (left), and feedback efficiency T] = P/C and the carbon dioxide emission intensity of energy 
c = E/a (right). Here, c represents the increase in atmospheric concentrations of CO2 E, per unit primary 
energy consumption a, that would be expected in a well-mixed atmosphere in the absence of terrestrial 
sink and source terms (1 ppmv CO2 = 2.13 Gt emitted carbon 132]). Dashed lines represent a least-squares 
first-order fit. Theoretical relationships between parameters are summarized in Table[T] 

Table 1 Mean observed and modeled quantities for the period 1970 to 2005. 



Parameter Functional dependence Observed mean Model mean 



energy efficiency growth 


d\nr\/dt 


1.05 %/yr 




carbonization growth 


d\nc/dt 


-0.29 %/yr 




feedback efficiency 


77 = )7,exp(^Mr) 


1.84 %/yr 




energy productivity growth 


dlnf/dt 


1.01 %/yr 


1.05 %/yr 


energy consumption growth 


dlna/dt = 77 


1.87 %/yr 


1.84 %/yr 


economic value growth 


dlnC/dt = rj 


1.82 %/yr 


1.84 %/yr 


economic production growth 


dlnP/dt = rj +d\nr]/dt 


2.88 %/yr 


2.89 %/yr 


CO2 emissions growth 


d In E/dt = rf+d\nc/dt 


1.56 %/yr 


1.55%/yr 
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